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Abstract. The supernova remnant S 147 harbors the pulsar PSR J 0538+2817 whose characteristic age is more 
than an order of magnitude greater than the kinematic age of the system (inferred from the angular offset of 
the pulsar from the geometric center of the supernova remnant and the pulsar proper motion) . To reconcile this 
discrepancy we propose that PSR J 0538+2817 could be the stellar remnant of the first supernova explosion in 
a massive binary system and therefore could be as old as its characteristic age. Our proposal implies that S 147 
is the diffuse remnant of the second supernova explosion (that disrupted the binary system) and that a much 
younger second neutron star (not necessarily manifesting itself as a radio pulsar) should be associated with S 147. 
We use the existing observational data on the system to suggest that the progenitor of the supernova that formed 
S 147 was a Wolf-Rayet star (so that the supernova explosion occurred within a wind bubble surrounded by a 
massive shell) and to constrain the parameters of the binary system. We also restrict the magnitude and direction 
of the kick velocity received by the young neutron star at birth and find that the kick vector should not strongly 
deviate from the orbital plane of the binary system. 
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star (NS)/SNR associations (e.g. Migliazzo et al. lMT^ . 

It is generally accepted that the supernova remnant 
(SNR) S 147 is associated with the radio pulsar PSR 
J 0538+2817. The only solid argument in support of this 
association is the positional coincidence of both objects: 

the pulsar is located (at least in projection) well within this paper, we propose an alternative explanation 

the extended shell of the SNR (Anderson et al. 1996). ^hc age discrepancy because PSR J 0538+2817 could 
The numerous estimates of the distance to the SNR are, the stellar remnant of the first supernova (SN) explo- 

in general, not inconsistent with the dispersion measure ^ion in a massive binary system and therefore could be as 

distance to the pulsar. The basic problem for the asso- oldas indicated by its characteristic age (cf. Morris et al. 

elation is the obvious discrepancy (Kramer et al. [Ml IM). Our proposal implies that S 147 is the diffuse rem- 

between the kinematic age of the system of - 3 x lO" yr ^ant of the second SN explosion (that disrupted the binary 

(estimated from the angular offset of the pulsar from the system) and that a much younger second NS (not neces- 

geometric center of the SNR and the pulsar proper mo- warily manifesting itself as a radio pulsar) should be asso- 

tion) and the characteristic (spin-down) age of the pulsar «ated with S 147. In Sect. 2 we review the existing obser- 

of ~ 6 X 105 yr. To reconcile these ages one can assume national data on the system PSR J 0538+2817/SNR S 147. 
that the pulsar was born with a spin period close to the Sect. 3 we suggest that the progenitor of the SN that 

present one (Kramer et al. IM see also Romani & Ng ^rmed S 147 was a Wolf-Rayet (WR) star and that the 

mm . This assumption is often exploited to explain the ^N explosion occurred within a wind bubble surrounded 

by a massive shell. In Sect. 4 we consider the possibility 

that PSR J 0538+2817 is the remnant of the first SN ex- 

* Address for correspondence: Krasin str. 19, ap. 81, Moscow plosion in a massive binary. Sect. 5 deals with some issues 

123056, Russia (vgvaram@sai.msu.ru) related to the content of the paper. 
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Fig. 1. The Hq image of the supernova remnant S 147 
(Drew et al. 120051 reproduced with permission of 
the IPHAS collaboration). Position of the pulsar PSR 
J 0538+2817 is indicated by a cross. The line drawn in 
the east-west direction shows the bilateral symmetry axis 
(see text for details). North is up, east at left. 

2. The SNR S 147 and PSR J 0538+2817: 
observational data 

2.1. General structure of S 147 

S147 (also G 180.0-1.7, Simeisl47, Shajn 147, etc) is a 
shell-type SNR with a diameter of ~ 3°. The optical im- 
age of S 147 presented by van den Bergh et al. I)1973|l shows 
a filamentary shell with a sharp circular boundary to the 
south. The north boundary of the shell is less regular: it 
consists of a long arc stretched in the east-west direction 
and two lobes protruding beyond the arc in the northeast 
and northwest directions for about one third of the charac- 
teristic radius of the SNR. There are some indications that 
the northern half of S 147 expands somewhat faster than 
the southern one (see Fig. 2 of Lozinskava ll976|l . The east 
and west edges of the SNR show signatures of blow-ups 
(more obvious to the east), that makes the SNR some- 
what elongated in the east-west direction. The blow-ups 
are seen more prominently in the recent excellent Hq im- 
age of S 147 presented by Drew et al. H2005I see Fig. EJ^ 
The 1.6 GHz image of S 147 by Kundu et al. 1(1^80) shows 
a break between the south and north halves of the SNR, 
typical of SNRs with bilateral symmetry. The bilateral 
axis of S 147, defined by the radio brightness distribution, 
is parallel to the east-west elongation of the SNR's shell. 

Optical and ultraviolet observations of S 147 (e.g. 
Lozinskava II 97KI Kirshner & Arnold . 1979;: Phillips et al. 
I1981|l suggest that the expansion velocity of its shell is 
~ 80 — 120kms~^, that implies that the SNR has already 
entered the final (momentum-conserving) stage of evolu- 
tion (see, however. Sect. 3). The same conclusion can be 
derived from the good positional agreement between sev- 
eral optical and radio filaments (Sofue et al. 119801 see also 
Fiirst & Reich [T^SSjl and from the non-detection of X-ray 
emission from the SNR's shell (Souvageot et al. I1990|l . 

2.2. PSR J 0538+2817 

PSR J0538-f2817 was discovered in the untargeted pul- 
sar survey conducted with the Arecibo radio telescope 
(Anderson et al. I199fi|l . which covered a 1° strip in the 
north half of S 147. The pulsar is located ~ 40' northwest 
of the geometric center of the SNR. The spin period P 
0.143 s and the period derivative P ~ 3.67 x 10~^^ss~^ of 
the pulsar (Anderson et al. I1996|l yield the characteristic 
age T — P/2P ~ 6.2 x 10^ yr, the magnetic field strength 
B = 3.2 X IQi^ {PP f/'^ ~ 7.3 X 10" G and the spin-down 

^ For other recent images of S 147 see 
:/ /www.skyfactory.org/simeisl47/simeisl47.htm| 



luminosity \E\ — Att'^IP ^P ~ 4.9 x lO^^'ergs ^, where 
/ ~ 10"*^ gcm^^ is the moment of inertia of the pulsar. 

The Chandra X-ray Observatory revealed a compact 
nebula surrounding PSR J 0538-^2817 (Romani & Ng 
12003(1 . Romani & Ng (,2003)) interpreted the nebula as an 
equatorial torus produced by the relativistic pulsar wind, 
similar to those observed around several other rotation- 
powered pulsars located in the high pressure interiors 
of their associated SNRs (e.g. Ng & Romani Ng 
& Romani H2004|l found that the symmetry axes of the 
toroidal nebulae (it is believed that they coincide with 
the pulsar spin axes) show a trend (most prominent in 
the case of the Crab and Vela pulsars) toward alignment 
with the proper motion vectors of the pulsars. Based on 
this trend, Romani & Ng l|20()3|l predicted the direction 
of the proper motion of PSR J 0538+2817 at a position 
angle of ~ 334° (measured north through east), which is 
consistent with the direction from the geometric center of 
S 147 to the present position of the pulsar (see, however. 
Sect. 5). 

Kramer et al. (I2003|l used timing observations of PSR 
J 0538+2817 to measure its proper motion in the eclip- 
tic coordinates, fi\ = —41 ± 3masyr~^ and fi/^ — 47 ± 
57masyr~^, with a median value /i = 67^22 ™^^sy^~^ ^^^'^ 

+28° 

a position angle of 311°^°o- The position angle is consis- 
tent with the movement of PSR J 0538+2817 away from 
the geometric center of S 147. Kramer et al. H2003|l con- 
sider this result as strong evidence that the pulsar is asso- 
ciated with the SNR (cf. Sect. 2.3). The large error in the 
proper motion measurement in the latitudinal direction, 
however, allows the possibility that the pulsar trajectory 
is significantly offset (up to ~ 0.5°) from the geometric 
center of S 147. 

2.3. Age of the system PSR J 0538+2817/SNR S 147 

For a long time it was believed that S 147 is one of the old- 
est evolved SNRs in the Galaxy. This belief was based on 
the age estimates of 1 — 2 x 10^ yr derived with help of the 
Sedov- Taylor solution or relationships for radiative blast 
waves (e.g. Sofue et al. lTMTI Kundu et al. 17801. It was 
further supported by the discovery of the associated pul- 
sar PSR J 0538+2817 of comparable (characteristic) age 
(Anderson et al. ll996|l . 

The subsequent studies of the pulsar, however, lead 
to a substantial reduction of the age of the system PSR 
J 0538+2817/SNR S 147. The proper motion measurement 
of PSR J 0538+2817 combined with its angular offset from 
the geometric center of S 147 yields the kinematic age, 
ikin, of only 3 x 10* yr (Kramer et al. I2003|l . Kramer 
et al. H2003|l pointed out that "since SNRs typically fade 
away after 100 000 yr, pulsars genuinely associated with 
SNRs are necessarily young" . To reconcile the discrepancy 
between tkin and r, they suggested that the pulsar was 
born with a spin period, Pq ~ 0.139 s, close to the present 
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one (cf. Romani & Ng I2(JU3P and therefore its true age 
is ^ T. They also admit the possibility that S 147 was 
formed by the SN explosion in a low-density bubble blown- 
up by the SN progenitor's wind (cf. Reich et al. i2003y) . i.e. 
the Sedov- Taylor solution cannot be used to estimate the 
parameters of the SN blast wave, and particularly the age 
of the SNR (see Sect. 3). 

We agree that S 147 could be the result of a cavity SN 
explosion (see Sect. 3) and note that in this case (i) the 
pulsar birthplace could be significantly offset from the ge- 
ometric center of the SNR due to the proper motion of 
the SN progenitor star, and therefore {ii) the proper mo- 
tion vector should not necessarily point away from this 
center (e.g. Gvaramadze 2002, Bock & Gvaramadze I2D021 
Gvaramadze 2004). Thus, (i) the true age of the sys- 
tem PSR J0538-t-2817/SNR S 147, t^ys, could be either 
< or > than i^in, and {ii) the association between PSR 
J 0538-1-2817 and S 147 could be genuine even if more ac- 
curate proper motion measurements will not prove that 
the pulsar is moving away from the geometric center of 
the SNR. 

Note also that the actual age of the system PSR 
J0538-f2817/SNR S147 is not fundamental for the fur- 
ther content of the paper (see, however. Sect. 5). The only 
thing we will assume in the following is that isys (or the 
age of S 147) is much smaller than r and that r is nearly 
equal to the true age of the pulsar. Even in the case of 
maximum possible offset of the SN blast center from the 
geometric center of S 147 and the current position of PSR 
J0538-f2817, Uys ^ 10^ yr. 

Now we discuss whether or not the spectral char- 
acteristics of the X-ray emission observed from PSR 
J 0538-f 2817 are at variance with our assumption that this 
pulsar is as old as indicated by its characteristic age. Thus 
one should check whether or not the effective temperature 
of the pulsar derived from model fits of its X-ray spectrum 
agree with the theoretical expectations. 

It is believed that the X-ray emission of middle-aged 
(~ 10^ yr) rotation-powered pulsars is predominantly of 
thermal origin and consists of two components. The first 
one, the hard thermal component, presumably originates 
from the hot polar caps heated due to the pulsar activity, 
and the second one, the soft thermal component, emerges 
from the rest of the surface of the pulsar. For older 10^ 
yr) pulsars, the stellar surface could be too cold to be 
observable in X-rays, so that the X-ray emission of these 
NSs comes mainly from the polar caps. 

The blackbody fits to the X-ray spectra of 
PSR J 0538-1-2817 obtained with the Chandra X-Ray 
Observatory (Romani & Ng l2()()3(l and the XMM-Newton 
(McGowan et al. I2003|l give the effective temperature 

^ Another possible way to reduce the age of the pulsar is 
to assume that the pulsar's magnetic field grows exponentially 
with a characteristic time-scale of ~ 4r/(n — fi) ~ 2 x 10** 
yr, where n = 3, n — vvjiP"^ and v,v and v are, respectively, 
the pulsar spin frequency and its two first time derivatives 
(measured by Kramer et al. I2()()3^ . 



Teft ~ 2.0 - 2.5 X 10^ K and the effective radius i?off ~ 1-3 
km of the spherical emitter. The inferred temperature is 
too high to be consistent with the temperatures predicted 
by the standard cooling models for NSs of age ^ t (e.g. 
Yakovlev & Pethick f2004|l . while small i?eff suggests that 
the X-ray emission is associated with regions much smaller 
than the entire surface of the NS. A possible interpretation 
of these fits is that PSR J 0538-1-2817 is a transition ob- 
ject between the middle-aged pulsars and the older ones, 
and that the observed X-ray emission is produced by the 
polar caps heated by the bombardment of relativistic par- 
ticles streaming down from the pulsar acceleration zones 
(McGowan et al. lM7S|l . 

Alternatively, the above temperature inconsistency 
could be considered as an indication that the surface 
of PSR J 0538-^2817 is not a perfect blackbody emit- 
ter and therefore its X-ray spectrum should be treated 
in the framework of a NS atmosphere model (with or 
without a mag netic field). McGowan et al. pMfl fitted 
the XMM-Newton spectrum of PSR J 0538-^2817 with 
a nonmagnetic hydrogen atmosphere model. They fixed 
the effective radius at values typical of NSs and found 
Teff ~ 0.6 - 0.7 X 10^ K and d ~ 0.3 - 0.4 k pc, wh ere d 
is the distance to the pulsar. Romani & Ng H2()03|l used 
the Chandra data to fit the spectrum of the pulsar with a 
magnetic hydrogen atmosphere model and derived about 
the same effective temperature, Teff ~ 0.65 x 10^ K; they 
also found i?efi = 13 km for d fixed at 1.2 kpc (the fig- 
ure based on the dispersion measure of PSR J 0538-1-2817; 
see Sect. 2.4). Although the temperatures derived in both 
models better agree with the standard cooling curves, 
they are still somewhat higher than the predicted ones. 
McGowan et al. (2003| suggested that their fit should be 
ruled out since it implies a much lower distance to the 
pulsar than the generally accepted value of 1.2 kpc (see, 
however. Sect. 2.4). 

Another possibility is that the pulsar PSR 
J 0538-1-2817 belongs to a class of very slowly cool- 
ing low-mass NSs with strong proton superfluidity in 
their cores, whose cooling curves lie above the basic 
standard cooling curve (see Yakovlev & Pethick i2004i and 
references therein). 

Thus we conclude that the existing X-ray data do 
not contradict the possibility that the true age of PSR 
J 05384-2817 is ~ r. 

2.4. Distance to the system PSR J 0538+2817 /SNR 
S147 

Most of the existing distance estimates for S 147 are based 
on the highly unreliable empirical relationships between 
surface brightness and linear diameter for SNRs. These 
estimates, ranging from d ~ 0.7 kpc (Milne I1970|l to 
~ 1.6 kpc (Sofue et al. 11980(1 . cannot be restricted from 
H I absorption measurements since S 147 lies towards the 
Galactic anticenter. A possible way to constrain the dis- 
tance to S 147 comes from the study of absorption lines 
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in spectra of stars located along the line-of-sight towards 
this extended SNR. Numerous observations (e.g. Phillips 
et aLUPl Philhps & GondhalekarTMl see also Sallmen 
& Welsh i200 4 ) have revealed the existence of high- velocity 
gas (associated with the SNR's shell) towards two stars, 
HD36 665 and HD37 318, located, respectively, at ~ 0.9 
and ~ 1.4 kpc (these figures were derived on the basis 
of spectral types, visual magnitudes and color excesses of 
the stars). Thus the distance to S 147 is ^ 0.9 kpc, that is 
consistent with the distance of ~ 0.8 kpc derived by Fesen 
et al. H1985fl from the interstellar reddening to the SNR. 
Note that the latter two estimates should be somewhat re- 
duced if the reddening to S 147 and the background stars 
is enhanced by the dust associated with the SNR's shell 
(cf. Gondhalekar & Phillips [1980|I . There are also several 
stars in the direction of S 147 with parallaxes measured 
with Hipparcos. The most distant of them, HD37367, is 
located at 0.36lg Jg kpc. The absorption spectrum of this 
star does not show high-velocity lines (Sallmen & Welsh 
I2()()4|l . that suggests that HD37 367 is a foreground star 
and puts a lower limit on the distance to S 147 of ~ 0.4 
kpc. 

An indirect estimate of the distance to S 147 comes 
from its association with PSR J 0538-1-2817. The disper- 
sion measure of the pulsar, DM~ 40pccm~'^ (Anderson 
et al. 1996), and the Cordes & Lazio ^U^^ model for 
the distribution of Galactic free electrons yield a distance 
to the system of ~ 1.2 ± 0.2 kpc, that is somewhat greater 
than the upper limit on the distance to S 147 of ~ 0.9 kpc. 
Note that the Cordes & Lazio model does not take into 
account a possible contribution to the dispersion of the 
pulsar's signal from the ionized material associated with 
the shell of S 147. If the excess dispersion measure due to 
the SNR's sheh ADM Z lOpccm-^ (see Sect. 3), than 
the distance to the pulsar is ^ 0.9 ± 0.2 kpc (Cordes & 
Lazio 120121) • 

In the following we allow d to vary in a wide interval 
from 0.4 to 0.9 kpc, and sometimes use the figure of 1.2 
kpc (accepted in the majority of recent papers devoted to 
PSR J 0538-1-2817 and S 147) to compare our results with 
those of other authors. The uncertainty in the distance, 
however, does not affect the main results of the paper. 



3. S 147 as the result of SN explosion within a 
WR bubble 

Let us now show that the small (kinematic) age of S 147 
and the low expansion velocity of its shell cannot be rec- 
onciled with each other if one assumes that the SN blast 
wave evolves in a homogeneous, uniform medium, i.e. if 
one describes the evolution of S 147 in the framework of 
the standard Sedov- Taylor model or models for radiative 
blast waves. 

According to the Sedov- Taylor model, the expan- 
sion velocity of the SN blast wave is given by vs-t = 
0.4i?sNR/isNR, where i?sNR and tsNR are, respectively, the 



See also ,http: / / rsd-www.nrl.navy.mil /7213 /lazio / ne_model/^ 



radius and the age of the SNR. The use of this model im- 
plies that the SN blast center coincides with the center of 
the SNR, and therefore tgNR = t^i^. For the angular radius 
of S 147 of ~ 1?5 and ^snr — 3 x 10'' yr, one has ws-t — 
340 di kms^^, where di is the distance to the SNR in units 
of 1 kpc. This estimate agrees with the observed expansion 
velocity ■ysNR ^ 80-120 kms"' if di = 0.24-0.35, that is 
for distances smaller than the lower limit given in Sect. 2.4. 
Moreover, the Sedov- Taylor solution implies the following 
estimate of the number density of the ambient interstel- 
lar medium: tijsm — 0.24 dj^^ cm^'^, i.e. ~ 300 cm^'^ for 
di — 0.24 or ~ 50 cm~'^ for di = 0.35. Both values of nisM 
are inconsistent with the fundamental assumption of the 
Sedov- Taylor model that the SN blast wave is adiabatic. 

The use of radiative models (e.g. Cioffi et al. 119881 
Blondin et al. I1998|l does not improve the situation. For 
example, using the relationships given in Cioffi et al. 
1)19881 see their Eqs. (3.32-3.33)), one has that vsnr could 
be consistent with the theoretical value if njsM — 10 cm~^ 
(for usNR — 120kms~^) or njsM — 30 cm^'^ (for usnr — 
80kms^^). In turn, these density estimates imply the dis- 
tance to the SNR of, respectively, ~ 0.41 and ~ 0.32 kpc. 

An alternative to these models is the possibility that 
S 147 is the remnant of a SN which exploded within a 
low-density bubble surrounded by a shell (created by the 
stellar wind of the SN progenitor). In this case, the ex- 
pansion velocity of the SN blast wave could be small even 
for young SNRs (see below). We believe that S 147 is the 
result of a SN explosion within the bubble blown-up dur- 
ing the WR phase of evolution of the SN progenitor star. 
Our belief is based on the interpretation of the general 
structure of S 147, whose shell is characterized by bilat- 
eral symmetry and is elongated along the symmetry axis. 

The bilateral and elongated appearance of some SNRs 
implies that the regular interstellar magnetic field is in- 
volved in shaping their shells. However, it was recognized 
long ago that the tension associated with the interstellar 
magnetic field cannot directly affect the shape of a typical 
SN blast wave to cause it to be elongated (e.g. Manchester 
11987(1 . On the other hand, the regular magnetic field could 
affect the symmetry of large-scale structures created in the 
interstellar medium by virtue of the ionizing emission and 
stellar wind of massive stars - the progenitors of most 
of SNe. For example, the elongated bilateral SNRs could 
originate if the SN blast waves in these SNRs take on the 
shape of the wind bubbles blown-up by the SN progenitor 
stars during the main-sequence phase and distorted by the 
surrounding regular magnetic field (Arnal '19921 Gaensler 
^1998). But the relatively small size of S 147 argues against 
a SN explosion within the main-sequence bubble and sug- 
gests that the pre-existing bubble was rather blown-up 
during the (much shorter) WR phase. Additional support 
for this suggestion comes from the low expansion veloc- 
ity of the SNR's shell, which could be naturally explained 
if the wind bubble was surrounded by a massive shell - 
the distinctive feature of WR bubbles (see Gvaramadze 
I 2004|l . In the presence of the regular interstellar magnetic 
field the structure of the WR shell is modified in such a 
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way that the density distribution over the shell acquires 
an axial symmetry with the minimum column density at 
the magnetic poles (see Gvaramadze |2UU4 and references 
therein). A good example of a bilateral WR shell is the 
nebula S308 around the WR star HD 50896 (see, e.g., 
Fig. Ic of van Buren & McCray '1988J. The subsequent 
interaction of the SN blast wave with the magnetized WR 
shell results in the origin of a bilateral SNR with two blow- 
ups along its symmetry axis. 

The low expansion velocity of the SNR's shell could be 
treated as an indication that the pre-existing (WR) shell 
was massive enough. Numerical simulations by Tenorio- 
Tagle et al. ([TMT|) showed that the SN blast wave merges 
with the wind-driven shell and the resulting SNR en- 
ters into the momentum-conserving stage (i.e. usnr ~ 
100 kms^^) if the mass of the shell is ^ 50 Mej, where Mcj 
is the mass of the SN ejecta. Assuming that the radius of 
the WR shell i?wR - 20 pc (rf = 0.9 kpc) and M^; ~ 4Mo 
(see Sect. 4), one has the number density of the ambient 
interstellar medium nisM > 0.2 cm^^. This estimate could 
be further constrained if one uses the result by Tenorio- 
Tagle et al. H1991|l that the reaccelerated wind-driven shell 
(now the shell of the SNR) acquires ^ 10% of the initial 
SN energy Eq = lO^i erg, i.e. vsnr ^ (0.2i;o/Af,hoii)'/', 
where Mshdi = (47r/3)i?^jiPisM, Pism = 1.4mHnisM, 
and mu is the mass of the hydrogen atom. For vsnr = 
lOOkms"^, one has nisM — LOcm^"^. The latter fig- 
ure could be used to estimate the excess dispersion mea- 
sure caused by the SNR's shell, ADM ~ i?sNR7iisM/3 — 
8(d/0.9kpc)^2p(,cjj^-3^ ^jjgj.g ^^^^ ^ 23 (d/0.9kpc) pc 

is the characteristic radius of the SNR (cf. Sect. 2.4). This 
excess could be much larger if a large-scale deformation 
of the SNR's shell (caused by the development of the 
Richtmaier-Meshkov and Rayleigh- Taylor instabilities in 
the reaccelerated wind-driven shell; see Gvaramadze J^£9 ) 
increases the line-of-sight extent of the ionized gas against 
the pulsar (cf. Gvaramadze I2()()rjl . 

The north-south asymmetry of S 147 could be under- 
stood if the WR shell was swept up from the medium with 
the density growing to the south (i.e. perpendicular to the 
orientation of the local regular magnetic field, implied by 
the bilateral symmetry of the SNR's shell). In this case, 
the SN blast wave merges with the south (more massive) 
half of the WR shell and this part of the SNR repeats the 
circular shape of the pre-existing shell. In the north direc- 
tion the blast wave merges only with the equatorial part 
of the WR shell (where the column density is enhanced 
due to the magnetic effect mentioned above) and overruns 
the less massive segments of the shell in the northeast and 
northwest directions, thereby producing two lobes in the 
north half of S 147. Our interpretation of the north-south 
asymmetry is testable. It would be interesting to measure 
(e.g. with a Fabry-Perot interferometer; see Lozinskaya 
I1976|l the expansion velocity of the northeast and north- 
west segments of S 147 to check whether or not it exceeds 
the characteristic expansion velocity of the SNR's shell of 
- 100 kms-^ 



Note that the above arguments in favour of cavity SN 
explosion should be relevant independent of whether S147 
is the result of the SN explosion of a single (WR) star or 
of the (second) SN explosion in a massive binary system. 

4. PSR J 0538+2817 as the remnant of the first 
SN explosion in a massive binary 

4.1. Constraints on the parameters of the binary 
system 

In Sect. 3 we suggested that the progenitor of the SN that 
created S 147 was a WR star, i.e. a massive star with 
the zero-age main-sequence (ZAMS) mass ^ 20 Mq (e.g. 
Vanbeveren et al. 11998(1 . Let us assume that this SN ex- 
plosion was the second one in a massive binary and that 
the pulsar PSR J 0538-^2817 is the remnant of the first SN 
explosion. The latter assumption implies that the ZAMS 
mass of the first SN progenitor was < 25 — 30 Mq (more 
massive progenitors produce black holes). We assume also 
that PSR J 0538-1-2817 is as old as indicated by its char- 
acteristic age, i.e. the SN explosions were separated by 
a time scale of ~ r. From this it follows that the ZAMS 
masses of the binary components were nearly equal to each 
other (cf. Bethe & Brown .19981 Vlemmings et al. IMlljl . 
so that it is likely that the second SN explosion also forms 
aNS. 

The spin characteristics (P and P) and the (inferred) 
magnetic field of PSR J 0538-1-2817 are typical of non- 
recycled pulsars. One can conclude therefore that the bi- 
nary system was sufficiently wide so that the stellar wind 
of the massive companion star did not appreciably affect 
the evolution of the pulsar, i.e. the standoff radius of the 
pulsar wind, r^, was larger than the accretion radius, ra 
(cf. Illarionov & Sunvaev ll985|l : 

r,{t) = i\E{t)\/A7:p^cviy/^ > = 2GMp/vl , (1) 

where c is the speed of light, pw = Mw/47ra^fw, M^, and 
are, respectively, the mass-loss rate and wind velocity 
of the companion star, a is the binary separation, G is the 
gravitational constant, and Mp = 1.4 Mq is the mass of 
the pulsar. From Eq. one has that a should be larger 
than some critical value given by the following relation- 
ship: 

2GMp ( M^v^c\ 

where |^*| = 32?:^ fi^ / 3c^ and = (P^ - 

IGtt^ fPtsys / ly are the spin-down luminosity and the 
spin period of the pulsar at the moment of the second SN 
explosion, /i = BR^ is the magnetic moment of the pulsar 
(we assume that fi — const), i? = 10 km is the radius of 
the pulsar, and isys ~ 3 x 10* yr. 

One can envisage two possible situations at the mo- 
ment of the first SN explosion. First, the companion star 
was a red supergiant and its extended envelope survived 
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the passage of the SN blast wave. Second, the convec- 
tive envelope of the companion red supergiant star was 
blown up by the SN ejecta to leave a bare He core (i.e. 
a WR star) or the companion star had already entered 
into the WR phase. In the first situation, = M^^*^ ~ 
lO~^M0yr~^ and = u^^^ ~ lOkms"^, so that one 
has from Eq. (0) that ~ 10^ Rq. In such a wide system 
the stripping and ablation of the red supergiant envelope 
can be neglected. In the second situation, Mw = ~ 
IQ-^Moyr-i and = u^R ~ 2 000kms-\ so that one 
has Ocr — 35i?0. Thus if a > Ocr, the pulsar was active 
during the whole period between the two SN explosions 
and its spin characteristics were not appreciably affected 
by the wind of the companion star; the latter implies that 
the true age of the pulsar is equal to t (provided that 
« P)- 

4.2. Origin of the pulsar peculiar velocity 

Proper motion measurement of PSR J 0538+2817 by 
Kramer et al. (2003) yields a pulsar transverse velocity 
of Vp^\ = 194 ± 14 di kms"""^, if one uses the proper motion 
in the ecliptic longitude only, or Vp = 318^^qq di kms^^ 
for the composite proper motion. For di = 0.4 — 0.9, one 
has Vp^x =i 80 - 170kms-i and Vp ~ 130 - 290 kms'^ 
Using the estimates of derived in Sect. 4.1, one can 
check whether or not Vp^x and Vp are consistent with the 
velocity, v^^, of the old NS released from orbit by the 
second SN explosion. 

In the case of a symmetric SN explosion and a circular 
binary orbit one has 

where m = M/Mp and M < 5 — 6Mq (e.g. Vanbeveren 
et al. 11998(1 is the mass of the pre-SN star. It is clear that 
^NS inconsistent with both Vp x and Vp if the second SN 
exploded after the red supergiant phase (i.e. a > 10^ Rq)- 
In the second situation (i.e. a > 35 Rq), one has from 
Eq. © that v^^ < 80 - llOkms^i. This estimate shows 
that f^'g is inconsistent with Vp, but could be equal to 
Vp_x if one adopts the smallest values of a and di and the 
largest one of M , and provided that the radial component 
of the pulsar velocity is small. 

One can, however, assume that the SN explosion was 
asymmetric so that the young NS received a kick veloc- 
ity, w, at birth. In this case, the new-born NS can impart 
some momentum to the old NS in the course of disinte- 
gration of the binary system (see Tauris & Takens'1998). 
The magnitude of the momentum depends on the angle, 
9, between the kick vector and the direction of motion of 
the exploding star, and the angle, 0, between the kick vec- 
tor and the orbital plane (see Fig. 1 of Tauris & Takens 
[TUMjl ^. An analysis of Eqs. (44)-(47) and (51)-(56) given 

o o o 

For = and < < 180 , the kick vector lies in the 
orbital plane and points to the half-sphere occupied by the old 
NS. 




e 



100 120 140 160 180 200 



Fig. 2. The dependence of the velocities of the old and 
new-born neutron stars (shown, respectively, by the solid 
and the short-dashed lines) on the angle between the kick 
vector and the direction of motion of the exploding star. 
The long-dashed line shows the best-fit transverse velocity 
at the upper distance 0.9 kpc. See text for details. 

in Tauris & Takens 1)1998(1 shows that the momentum im- 
parted to the old NS is maximum if 

0^9^,= arccos(— w/w) , (4) 

where v = [G{M + Mp)/a]^/^ is the relative orbital veloc- 
ity, and provided that the vector of the kick velocity does 
not strongly deviate from the orbital plane of the binary 
system; i.e. than the kick received by the second-born NS 
is directed almost towards the old NS. Fig. |21 illustrates 
how the direction of the kick affects the velocities of the 
old and new-born NSs released from the disrupted binary. 

o 

One can see that for ~ 6^ — 118 (we assume that 
w = 400kms^\A/ = SMq and a = 35 i?©, and that 

o , . 

(/) = ) ti^g is maximum (~ w), while wj^g" drops to a 
minimum value of ~ 0.25w. It is also seen that for 9 rang- 

o ° IH 

ing from ~ 115 to 2± 135 w^g is larger than the best-fit 
transverse velocity at the upper distance 0.9 kpc. Note 
that the maximum value of is almost independent of 
a (or v): u^'g'™'^'^ tends to w with increasing a^, while 
drops to 0. Fig. 121 shows how w^'g depends on 0; it is 
seen that w^'g'"^'''^ ^ w™'*'' if ^ 10 . The above consid- 
erations show that u^'g could be consistent with Up, but 
the binary separation must be as small as possible while 
avoiding pulsar recycling and the kick direction must be 
carefully tuned. For di — 0.9 and a ^ 35Rq the proba- 
bility of a favourable kick orientation is ^ lO^'^ (for the 
isotropic kick distribution) or ~ 10~^ (if kicks produced 
by SN explosions in binary systems are restricted close to 
the orbital plane; see Sect. 5). Note also that the smaller 
di the wider the range of angles for which > Vp and 
the larger the probability that the second-born NS will 
receive an appropriately oriented kick. 

^ Vlemmings et al. 12UU4II draw an erroneous conclusion that 
in a wide binary disrupted after the second (asymmetric) SN 
explosion the old NS could be accelerated to a velocity ~ w, 
while the second-born NS to a velocity several times exceeding 
w. 
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Fig. 3. The dependence of the maximum value of the ve- 
locity of the old neutron star on the angle between the 
kick vector and the orbital plane. 




Fig. 4. The angle between the velocity vector of the old 
neutron star and the orbital plane as a function of the 
angle between the kick vector and the orbital plane. 



5. Discussion 

In Sect. 4.2 we considered two possible explanations for 
the origin of the peculiar velocity of PSR J 0538+2817 
based on the idea that the pulsar could be the remnant 
of the first SN explosion in a massive binary. Proceeding 
from this we found that although the pulsar motion in 
the ecliptic longitude taken separately (cf. Hobbs et al. 
120051 Lewandowski et al. I2()()4|l could be consistent with 
the possibility that the binary was disrupted by a sym- 
metric SN explosion, the large composite proper motion 
measured by Kramer et al. (?; if confirmed by more accu- 
rate measurements) makes the case of an asymmetric SN 
explosion more plausible. In this case, we expect that the 
kick velocity received by the second-born NS is restricted 

o 

close to the orbital plane (i.e. ^ 10 ; cf. Wex et al. 
I2000II . Fig. ^ shows the angle, x, between v^^ [for 9 given 
by Eq. 10}] and the pre-SN orbital plane as a function of 

o o o o 

(p: X grows from (for = ) to ~ 40 (for ~ 10 ) 

o o 

and than gradually decrease to ~ 10 for (p = 90 . Thus 
if our explanation of the age discrepancy is correct, we 
do not expect any alignment between the pulsar spin axis 
and proper motion, although one cannot exclude that they 
are perpendicular to each other if the second SN explo- 



sion was symmetric about the orbital plane (i.e. ^ ~ ; 
cf. Cordes & Wasserman lTOMI Colpi & Wasserman 
Vlcmmings et al. 2004) and provided that the spin axis 
of PSR J 0538-1-2817 was perpendicular to this plane. To 
check whether this possibility was realized one must know 
the precise position angles of the pulsar proper motion 
and the symmetry axis of the X-ray nebula surrounding 
the pulsar. However, the large error in latitudinal direc- 
tion makes the position angle of the pulsar proper motion 
very uncertain (see Sect. 2.2), while the low photon statis- 
tics of the Chandra data does not allow us to infer the 
symmetry axis of the pulsar wind nebula unambiguously. 
We will discuss this point in more detail. 

The existing Chandra data on PSR J 0538-f2817 show 
that the nebula around the pulsar consists of a core and a 
dim halo (Romani & Ng l2003j) . These substructures could 
be fitted with ellipses with the semi-major axes of 2'.'5 and 

n o o 

9 , at the position angles of ~ 130 and ~ 60 . Romani 
& Ng (j2()03|l suggested that the halo "has a scale compa- 
rable to the expected wind shock radius for a confinement 
pressure appropriate to the interior of S 147" and used 
the geometry of the halo to infer the pulsar spin axis at a 
position angle of 334? ± 5?5 (that is consistent with the 
position angle of the pulsar composite proper motion). 
Their calculations of the shock radius were carried out in 
the framework of the Sedov- Taylor model and under the 
assumption of a spherically-symmetric pulsar wind. 

The angular radius of the toroidal shock can be ex- 
pressed as 



e ~ (|£:|/47racd2pi„)i/2 



(5) 



where a < 1 is a factor characterizing the anisotropy of 
the energy flux of the pulsar wind (for a << 1 most of 
the pulsar wind is confined to the equatorial plane, while 
a — I corresponds to the spherically-symmetric wind) , 

(7 - l)£^th 



Pi 



V 



(6) 



is the mean gas pressure inside the SNR, 7 is the specific 
heat ratio, Eth = SEq is the thermal energy of the SN 
blast wave, 6 < 1, and V — (47r/3)i?gi,^j^ is the volume 
of the SNR. For adiabatic gas (7 = 5/3) and an angular 
radius of the SNR of ~ 1?5, one has from Eqs. (O and l|njl 
that 



e ~ l'.'4(a<5)-l/2rf;/^ 



(7) 



For a Sedov- Taylor blast wave (i.e. S ~ 0.72) and as- 
suming a = 1 and di = 1.2 (Romani & Ng'2003'), one has 
Q = 1'.'8, that is consistent with the angular radius of the 
core of the nebula rather than with the scale of the more 
extended halo. Moreover, in Sect. 3 we showed that the 
Sedov- Taylor model cannot reconcile the small expansion 
velocity of S 147 with its low kinematic age and suggested 
that this SNR is the result of a SN explosion within a 
pre-existing bubble surrounded by a massive wind-driven 
shell. In this case, the evolution of the SN blast wave com- 
pletely differs from that of the Sedov- Taylor one. The pres- 
ence of the massive shell strongly increases the radiative 
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cooling of the SN blast wave, that results in a rapid de- 
crease of the thermal content of the SNR as compared with 
the Sedov- Taylor one. Numerical simulations by Tenorio- 
Tagle et al. ijTMTl see their Fig. 3) show that S ~ 0.4 
at tsNR = 3 X 10^ yr and further decreases to ~ 0.1 at 
^SNR = 10^ yr. Thus assuming that ^snr = ikin = 3 x 10* 
yr, one has from Eq. JT)) that Q < 2".l if the distance to 
S 147 is < 0.9 kpc, or ~ 2'.'4 if di = 1.2. These estimates 
also suggest that the toroidal shock could be associated 
with the core of the pulsar wind nebula. However, in the 
case of a cavity SN explosion tkin could be as large as ~ 10^ 
yr (see Sect. 2.3); in this case ,5 ~ 0.1 and 6 ~ 3'.'6 - 4'.'7 
(for di ranging from 0.9 to 1.2), that better agrees with 
the scale of the halo (cf. Romani & Ng 12003(1 . The agree- 
ment could also be achieved when the anisotropy of the 
pulsar wind is taken into account. Assuming that a = 0.1, 
one has from Eq. Q) that 6 < 6'.'6 and 7'.'6 for di = 0.9 
and 1.2, respectively. 

If the core indeed corresponds to the equatorial torus of 
the pulsar wind nebula (this could be checked with forth- 
coming deep Chandra observations^ or high-resolution ra- 
dio imaging) , than the position angle of the symmetry axis 
of the torus is exactly perpendicular to the direction of the 
pulsar composite proper motion (intcrferometric measure- 
ments of PSR J 0538-^2817 would be highly valuable for 
further restricting the position angle of the pulsar proper 
motion), i.e. just what is expected if the pulsar obtained 
its peculiar velocity by the disintegration of the binary 
system (with aligned angular momenta) by symmetric SN 
explosion or by SN explosion symmetric about the orbital 
plane (see above). Note, however, that the same spin-kick 
orientation could be produced by the SN explosion of a 
solitary massive star if the spin and the peculiar velocity 
of the pulsar are due to a single off-center kick (Spruit & 
Phinnev lT^Mjl . 

Now we discuss the problem of the second (young) NS 
possibly associated with S 147. The detection of such a NS 
would allow us not only to distinguish the case of a SN 
explosion in a binary system from that of a solitary SN 
explosion, but also would lend strong support to our ex- 
planation of the age discrepancy and would have a strong 
impact on our understanding of the kick physics. 

In the case of a symmetric SN explosion the angle, -0, 
between the velocity vectors of the old and new-born NSs 
depends only on the mass of the SN ejecta and is given by 
(Gott et al. lWTni see also Iben & Tutukov 11996,1 



arccos 



2Mr 



-1/2 



(8) 



where Mej = M - Mp. It follows from Eq. © that ip is 

o o 

always smaller than 90 ; e.g. -0 ~ 60 for the parameters 
adopted above. In the case of an asymmetric SN explosion 

o _ 

ip could be larger than 90 . Fig.Oshows how i/j depends on 
the direction of the kick vector (we assume here that — 

o 

). It is seen that ip is very sensitive to 0. For ranging 
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Fig. 5. The angle between the velocity vectors of the old 
and new-born neutron stars as a function of the angle 
between the kick vector and the direction of motion of the 
exploding star. 



from 115 to 135 (i.e. when v' 



old 
NS 



for 



^ 9„ [see Eq. (0}] and grows to ~ 60 —70 at the bounds 
of the range. Thus we expect that the young NS should 

o 

be located in a cone with a half-opening angle ^ 70 with 
the cone axis oriented along the proper motion of PSR 
J 0538+2817. We also expect that the young NS should 
be situated closer to the vertex of the cone (i.e. to the SN 



blast center) than PSR J 0538+2817 since v- 



now 
NS 



< V- 



,old 
NS 



(the 



same is true for the case of a symmetric SN explosion). 
Remind that the SN blast center could be significantly 
offset from the geometric center of S 147 (see Sect. 2.3). 

The most optimistic supposition is that the young NS 
is an ordinary (rotation-powered) pulsar with a favourably 
oriented radio beam. The youth of the pulsar implies that 
it should be more energetic than PSR J 0538+2817, and 
therefore could be easily detected somewhere to the south 
of the 1 strip covered by the Arecibo survey (most likely 
in the west half of the SNR). Another possibility is that 
the young NS is an off-beam radio pulsar or that it belongs 
to a class of radio-quiet NSs. In both cases one can expect 
that the young NS should be a sufficiently bright soft X- 
ray source to be detected with the i? T" All-Sky Survey 
(RASS). Indeed, the RASS Faint Source Catalog'' shows 
five sources in the proper place. The count rates of these 
sources (ranging from ~ 0.014 to ~ 0.039 cts^^) should be 
compared with the count rate in the ROSAT pass band 
expected for a NS of age of ~ 3 x 10** yr. Let us assume 
that the new-born NS is a spherical blackbody emitter of 
radius of 10 km and temperature of ~ 10^ K (predicted 
by standard cooling models). Then assuming an interstel- 
lar absorption column density of ~ 2.5 — 3 x lO^'^ cm~^ 
(e.g. McGowan et al. 12003 Romani & Ng lM)^ and us- 
ing PIMMS^ one has ~ 0.017 - 0.022 ROSAT PSPC 
ct s~^, that is, the figures comparable with the above count 
rates. Thus one cannot exclude that one of the five RASS 
sources is a young NS associated with S 147. If this NS is 
a rotation-powered pulsar, its position should be marked 



http:/ /www.xray.mpe. mpg.de/rosat/survey/rass-fsc/ 
See ,http://heasarc.gsfc. nasa.gov/cgi-bin/W3Browse/w3table. pH http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 
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by a pulsar wind nebula with a characteristic scale of at 
least an order of magnitude larger than that of the neb- 
ula around PSR J 0538+2817. Deep radio or X-ray obser- 
vations would allow us to verify the existence of such a 
nebula. 

Note that the RASS sources are about 1.5-4 times 
fainter than the pulsar PSR J 0538-1-2817 (also detected by 
the RASS; e.g. Sun et al. ll996|l . while the opposite situa- 
tion is expected if PSR J 0538-1-2817 is as old as indicated 
by its spin-down age and if its cooling follows the stan- 
dard cooling curves. The contradiction, however, could be 
removed if (as discussed in Sect. 2.3) the pulsar belongs to 
a class of slowly cooling NSs. 

One cannot exclude that the young stellar remnant 
collapsed into a black hole. In this case, the chances of 
finding this object are small. 
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